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The electric conductivity which reveals the zero gap semiconducting (ZGS) state has been 
investigated as the function of temperature T and life time r in order to understand the ZGS 
state in quarter-filled a-(BEDT-TTF)2l3 salt with four sites in the unit cell. By treating r as 
a parameter and making use of the one-loop approximation, it is found that the conductivity 
is proportional to T and r for fcsT 3> H/t and independent of T and r for ksT <C h/r. 
Further the conductivity being independent of T in the ZGS state is examined in terms of 
Born approximation for the impurity scattering. 
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In quasi-two-dimensional organic conductors, 
bis(ethylene dithiolo)tctrathiofulvalene (abbreviated 
as BEDT-TTF) salts, several electronic states, such as 
Mott inslator, charge ordering and superconductor, 1 ' 2 
have been obtained by varying temperature, hydrostatic 
pressure and uniaxial strain. Especially, the unconven- 
tional behavior with decreasing temperature has been 
found in the a-(ET) 2 I 3 salt (a phase of BEDT-TTF 
tri- iodide), which consists of four molecules per unit cell 
and 3/4- filled band (quarter filling of hole) with eight 
transfer energies. 3 For ambient pressure, a-(ET)2l3 
shows a metal-inslator transition at Tmi = 135 K 
followed by the charge ordering. 4 Under hydrostatic 
pressure, Tmi decreases and vanishes where, at 20 kbar, 5 
the resistivity becomes almost temperature-independent 
from 300 K to 1.5 K. The carrier density decreases from 
10 21 cm" 3 (at 300 K) to 10 15 cm" 3 (at 1 K) while 
the carrier mobility is interpreted to move from 10 _1 
cm 2 /(V-s) (at 300 K) to 10 5 cm 2 /(V-s) (at 1 K). It has 
been maintained that such a temperature-dependence 
of the carrier indicates a narrow gap semiconducting 
(NGS) state with the gap E g ~ lmeV. 5 The NGS state 
is also obtained at 10 kbar, 6 under the uniaxial strain 7 
along a-direction. 

The band calculation of a-(ET)2l3 shows a zero gap 
semiconducting (ZGS) state 8,9 instead of the NGS state 
when the uniaxial pressure along a- axis, p a , becomes 
large. The Fermi surface of a-(ET)2l3 salt at ambient 
pressure shows a semi-metallic state with the hole and 
the electron pocket, 5 and the band spectrum has the de- 
generacy at the bottom of the electron pocket with an 
incommensurate wave number ±fe . 10 Since k° is differ- 
ent from the location corresponding to the symmetry of 
the Hamiltonian, such a degeneracy is called as acciden- 
tal one according to Herring. 11 With increasing p a , the 
hole and electron pockets are reduced 12 and diminish, for 
p a > 3 kbar, with the resultant band dispersion, showing 
Dirac cones centered on ±fc . 10 Such a band dispersion 
and the corresponding density of states per spin and per 



volume are given respectively by 

£fc = ±hv\k - k°\ 

D(e) 



(1) 

close to the Fermi energy Sp, where v is the average of 
the velocity and e is the energy measured from ep. 

In terms of eq. (1), the carrier density n of the ZGS 
state at low temperatures is calculated as 

r 2 (fc B T) 2 



n = 2 



D(e)f(s)ds = 



1 



Trh 2 v 2 



(2) 



where /(e) = (e £ / fesT + l) -1 and T is the temperature. 
Equation (2) is consistent with that of the experiment. 5 ' 6 
However, it is unclear if the increase of the carrier den- 
sity is canceled by the decrease of the carrier mobility 
when the conductivity keeps constant with increasing 
temperature. In the present paper, the conductivity as 
the funcion of T and r is examined for the ZGS state 
of a-(ET)2l3- The calculation is performed by using the 
Green's function without the vertex corrections where 
the life time r in the self energy is treated as a parame- 
ter. 

In the system with multi-sites, we consider explicitly 
H for the transfer energy which is crucial to obtain 
the ZGS state, but treat implicitly H int for the inter- 
action which gives rise to the life time. The coupling 
to the external field, H ext is also added to calculate 
the conductivity. The total Hamiltonian is written as 
H = H + H int + H exU where 



H = 



n.n. 



Hext = -e^r-i • Ea\ a a ia = - A, 



(3) 



(4) 



and i,j(= 1, ••• ,Nl) denotes indices for the unit cell, 
forming a square lattice. The indices a, [3 denote the site 
in the unit cell. The quantity ti a: jp is the transfer energy 
for the nearest neighbor (n.n.) sites between (i, a) and 
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(j, (3), and a,i a is the annihilation operator of the conduc- 
tion electron at (i,a) site. In eq. (4), n is the position 
vector of the unit cell coupled with the electric field E 
where e is the electric charge. Since the equation of mo- 
tion for A is written as A — (i/H)[H, A] — J ■ E, the 
current density J p (/i being the direction of J) is given 
as 

e ^ ds al3 (k) t e ^ t 

= a E — a L«fe/s = ^ E < 7 '(fe)4 7 c fe7' 5 

(5) 

where 

^H^EW^"', (6) 

71. 71. 

<y(fe)= E d* ai {k)^^d Pl ,{k\ (7) 

af3=l ^ 

and V" is the area of a system. The directions x and y 
correspond to a and 6 in the Mori's notation, respectively 
(inset of Fig. I). 3 The quantities ak a and Cfc 7 denote the 
Fourier transform, a^a = V~ x l 2 X^=fi ai a e~ lk " ri , Cfc 7 = 
ELi^W^c" where the band index 7(= 1, • ■ ■ ,4). 
The quantity d ai (k) is the component of the eigenvector 
obtained by 

4 

^2(e af3 {k) - e F )dpj(k) = £fc 7 d a7 (fc), (8) 
/3=i 

where £fc 7 (£fci > £fc 2 > £fc3 > 6c4) is the eigenvalue 
of the multi-band with the state (fe,7). The matrix ele- 
ment £a>p(k) is given by eq. (6), which depends on eight 
transfer energies, tA, as shown in the inset of Fig. 1. 
The indices A = cl, c2, c3, c4,pl,p2,p3,p4 correspond to 
a2, a3, al, ol, 62, 61, 64, 63 , respectively, in the definition 
by Mori et al. 3 Since the ZGS state is 3/4-filling, the 
band of £fci is empty and the others are filled at T = 0. 

In terms of the linear response theory, the electric con- 
ductivity a**", (/i, v = x, y) is expressed as, 

a*"' = i[(Q* (w) - Q* (0))M^o • (9) 

The retarded Green's function Q^ u (uj) is obtained from 
the analytical continuation of Q^(iuji) where Q^(iuii) 
is the two-body Green's function (Tj is the ordering op- 
erator) of imaginary time U, given by 

1 f (i 

Q llv (m) = -y J (TiJ^UWe^dU, (10) 

where uii = l-KksTl with I being an integer. 

Equation (10) is calculated by adopting one-loop ap- 
proximation, i.e. without vertex correction. We treat the 
effect of Hint as a life time, r, in the Green's function, 13 

rP 

G 7 y(fc,ze m ) = -J (T 4 c fc7 (^)4 7 '(0))e 4£mtl d^ 

- (i£ m - Cfe 7 - E(ie m )) _1 5 77 < , (11) 

E(ze m ) = -ift(s g n(e m )/2T) ) (12) 

where e m = (2m + l)7rfcsT with m being an integer. 
The life time, r, is naively taken as a parameter where r 
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may come from the processes such as the scattering by 
clcctron-phonon interaction, impurity and fluctuation. 
Within the above approximation, eq. (10) is expressed 
as 

Q^i) = e ^f E < y (fcK' 7 (fc) 

x Gyy(k,i(e m + uji))G 71 (k,ie m ) 

p 2 . r°° 

= ^E<7'( fe K%( fe )/ 

X [i + iwj - f fc y + (ift/2-r) (e - Cfe 7 ) 2 + (ft/2r) 2 
1 (H/2t) 

+ e-iui- & 7 - {ih/2 T ) (e - £ fey ) 2 + (V2r) 2 

(13) 

After the analytical continuation (i.e., iw; — ► tko + iO) , 
the real part of (/i = v = x, y) per spin is calculated 

as, 

^ = ReK' t )=^< 7 „ 

77' 

(fi/2r) (fi/2r) 

X (e - ^ 7 ') 2 + (^/2r) 2 (e - £ fe7 ) 2 + (ft/2r) 2 ' 

(14) 

Based on the experimental data, we use the interpola- 
tion formula for the transfer energy, where tA (eV) at p a 
(kbar) is given by , 8, 12 

t A (Pa) =t A (l+K APa ). (15) 

The units of the energy and the pressure are taken as 
eV and kbar, respectively. In eq. (15), t c \ = 0.048, t c2 — 
-0.020, t c3 = -0.028, t ci = -0.028, t pl = 0.140, t p2 = 
0.123, t p3 = -0.025, t p4 = -0.062, K cl = 0.167, K c2 = 
-0.025, K c3 = 0.089, K c4 = 0.089, K pl = 0.011, K p2 = 
0.000, K p3 = 0.000 and if p4 = 0.032. The ZGS state is 
obtained for p a > 3. 10 We calculate a as the function of 
T and r where the area of an unit cell V/Nl, the electric 
charge e, Plank constant h and Boltzmann constant ks 
are set unity. 

Figure 1 depicts the temperature dependence of a for 
1/t = 0.001 and p a = 10 (ZGS state). The difference 
between a x and a v comes from the anisotropy of t A , 
where a x is larger (smaller) than a v for T > 0.14 (T < 
0.14) in the present choice of parameters. In Fig. 1, one 
finds cr^ cx T at low temperatures in contrast to the 
conventional conductivity given by the Drude's formula 
(i.e. a — ne 2 T /m with m being the electron mass). Thus 
it is incorrect that <r M in the ZGS state is propotional to 
the carrier number with n(oc T 2 ). The decrease of ct m at 
low temperatures comes from the density of states D{e), 
which is given by eq. (1) near the Fermi energy. 

Figure 2 shows the temperature dependence of a y for 
several choices of life time 1/t. In the high temperature 
region, the conductivity increases monotonically up to 
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Fig. 3. Main figure: temperature dependence of <r v (dotted line), 
Fig. 1. Temperature dependence of <j* and uV for ZGS state at ^ (soUd line)> ^^^led line) and a\ 2 (= c v 21 ) (dot-dashed 

p a = 10 (kbar), where 1/r = 0.001 (eV). The inset shows the lme) at pa = 1Q and 1/r = rj.OOl. Inset: the corresponding quan- 
structurc of Q-(ET) 2 I 3 on the conducting plane. titiy as thc function of 1/r at Pa = 10 and T = 00L 




10" 2 T[eV] 10" 1 



Fig. 2. Temperature dependence of cr y at p a = 10 where the in- 
set shows the 1/r dependence of a v at the temperature T = 
0.01(i),0.003(ii),0.001(iii), and 0.0003(iv) respectively. 



T = 0.1, and is large for small 1/r. In the low temper- 
ature region, a y is independent of the temperature This 
region is wider when 1/r is larger. In the inset thc 1/r 
dependence of a y is shown for T = 0.01, 0.003, 0.001 and 
0.0003. Thc a y is flat for large 1/r and the gradient of a v 
is about 1 for small 1 jr indicating that a y is propotional 
to r. The above result is different from the conventional 
metallic state in which the conductivity is always given 
by a — ne 2 T jra. The region in which the conductivity is 
independent of the life time is wide when the tempera- 
ture is low. 



In Fig. 3, the temperature dependence of the compo- 
nents of cr^ 7 , with (7,7') = (1,1), (1,2), (2, 2) and a y are 
shown for p a = 10 and 1/r = 0.001. The correspond- 
ing conductivity as the function of 1/r (T = 0.001) is 
also depicted in the inset. The diagonal components of 
a\ x and ct| 2 mainly contribute to a y , whereas the off- 
diagonal component crf 2 (— a 2i) contributes only at thc 
low temperatures. The fact that a\ 2 decreases monoton- 
ically in contrast to a\ x and cr| 2 m the region of high 
temperature is understood as follows. The carriers ex- 
cited thermally can participate in g\ x and er| 2 while only 
carriers close to Fermi surface give rise to a\ 2 and a 21 . 
For the 1/r dependence, a\ x and a 22 are propotional to r, 
and a\ 2 decreases for small 1/r. Such a behavior of o\ 2 is 
also understood by noting that the corresponding height 
in the r.h.s. of cq.(14) at e = becomes large for large 
1/T and large 1/r. Thus a\ 2 can contribute in the dirty 
case of Tr <g; 1. Compared with these four components 
of of-L, a 22 , af 2 , and a 22 , the other twelve components of 
cr^ 7 , are negligibly small. 

The behavior of the conductivity in Fig. 2 can be ex- 
plained from thc property of the density of states in eq. 
(14). Since D(lo) = J^k 6 ^ ~ = (M/ 2 ™ 2 ) with 
±£fc = V\Jk 2 + k% (v is an averaged velocity), eq. (14) 

with only the diagonal component of (7,7') = (1, 1) and 
(2, 2) is calculated as (x — 0e), 



a ~ — / du>deD{ijj) 



dm 

de 



(l/2r) 



(l/2r) 



1 

4^2 



dx 



1 + (2Tr|x|)tan- 1 (2Tr|a;| 



(16) 



cosh x + 1 

When T > 1/r, one can replace tan" 1 (27Y|x|) by tt/2 
and estimate eq. (16) as 

a ~ (log 2/7r)Tr ~ 0.220Tr. (17) 
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Fig. 4. Conductivities by Born approximation (dashed line) , 
which arc compared with those of Fig. 2 (solid line). 



Thus a is propotional to the temperature and the life 
time for T> 1/r. The results of eq. (17) is also derived 
by using Boltzmann's equation where the conductivity, 
<tb, is expressed as 




Actually eq. (16) becomes equal to eq. (18) for |e| » 1/r, 
indicating that the case T ^> 1/r corresponds to the 
classical limit. When T <C 1/r, eq. (16) is calculated as 

a ~ 1/(2tt 2 ) = 0.051, (19) 

which shows a being indepnedent of both temperature 
and life time. In both cases of T ^> 1/r and T <C 1/r, 
the conductivity does not depend on the velocity v. These 
two regions of a ~ const, and <r ~ Tr are seen in Fig. 2. 
The boundary between two the region of a ~ const, and 
that of ct ~ Tr is obtained as 1/r ~ 0.7T at p a = 10 and 
1/r ~ 0.8T at p a = 6. Equaions (17) and (19) exhibit 
the universal coefficient of Tr and the universal constant, 
respectively. These values multiplied by 2 are compared 
with those of Fig. 2, because eqs. (17) and (19) are calcu- 
lated in a single Dirac cone. For 0.001 < T < 0.005, the 
coefficient of Tr in Fig. 2 corresponds well to the uni- 
versal coefficient. At low temperature, however, the con- 
ductivity in Fig. 2 is somewhat larger than the universal 
constant, because the former contains the contribution 
from the off-diagonal components, (7, 7') = (1, 2), (2, 1). 

Now we discuss the relevance of the ZGS state to the 
conductivity of a-(ET) 2 i3 , which is almost independent 
of temperature. In terms of the mobility and the carrier, 6 
the bulk conductivity is estimated as a ~ 3.2 x 10 2 O — 1 
cm -1 . Then one obtains 1/r ~ IK at room temparature 
by using naively the formula a = ne 2 r/m with the bare 
mass m. Since the present result of the temperature de- 
pendence for T > 10K (i.e., Tr » 1 ) is not consistent 
with that of the experiment, we calculate the life time 
of the impurity scattering by taking account of the ZGS 
state. In terms of the self-consistent Born approxima- 



tion, 13 the self-energy in eq. (12) is replaced by 

E(fe m ) = -ihe m /(2 T W D ), (20) 

for \e m \ < Wd, where Wd (~ 0.01 at p a — 10 kbar) de- 
notes the band width of the Dirac cone. Note that the 
frequency dependence similar to eq. (20) has been ob- 
tained for the calculation in the graphite at T = 0. 14,15 
Thus we obtain the conductivity in which the tempera- 
ture dependence is reduced as shown in Fig. 4 indicating 
a tendency being consistent with the experiment. How- 
ever there remains a problem of the magnitude of a which 
comes from the stacking of the two-dimensional sheet of 
the ET salt. From the balk conductivity and the lattice 
constant of c-axis (as 1.7 x 10~ 7 cm), 12 the conductivity 
per one sheet and spin, is estimated as 0.11 in the unit 
of e 2 /h(= 2.44 x 10~ 5 fi _1 ), which is utilized for the con- 
ductivity in the present calculation. This magnitude is 
is much smaller than the numerical result of Fig. 4 al- 
though the origin of such a discrepancy is not clear for 
the moment. 

Here, we show the validity of the extrapolation formula 
cq. (15) even for 10 kbar by comparing the coeffients of 
the density of states in cq. (1). The experimental value of 
the carrier density at T is given by n ~ 10 15 T 2 cm -3 and 
then the density per sheet is given by n ~ 1.7 x 10 8 T 2 
cm~ 2 . By comparing the coefficient with that of eq. (2) 
( lattice constants are a = 9.2 x 10 -8 cm and b — 10.8 x 
10~ 8 cm ), 12 the coefficients of e in eq. (1) is estimated as 
8 x 10 cV~ 2 , which corresponds well to that the present 
calculation ~ 6 x 10 cV~ 2 . 

Finally we note that the three dimensional effect on 
the ZGS state. Since the small interplane hopping is ex- 
pected in the a-(ET) 2 i3 salt from the first principle cal- 
culation, 16 the metallic state with small Fermi surface 
appears. 17 Then the property of the ZGS state could be 
found at temperature larger than the interplane hopping. 
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